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The multichannel quantum defect theory (MQDT) can be reinterpreted as a quan- 
tum Poincare map in representation of angular momentum. This has two important 
Q^ . implications: we have a paradigm of a true quantum Poincare map without semi- 

classical input and we get an entirely new insight into the significance of MQDT. 
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In recent years there has been a rapidly growing interest in the quantum 
Poincare map (QPM) [1-10], i.e. the quantization of a classical Poincare map, 
for a time independent Hamiltonian system. Bogomolnyfl] started out with a 
semi-classical formulation. Among other things he shows that unitarity of the 
representation is reached in the limit h — > 0, while this is only approximate 
for finite H [8]. Prosen [4] gives an elegant general solution to the unitarity 
problem at the expense of obtaining an infinite matrix for the QPM. The semi- 
classical approach common to most discussions causes a number of problems 
that make the use of this new and powerful tool a little obscure. With other 
words, the quantum Poincare section implicitly defined by Bogomolny, lacks a 
paradigmatic example where a quantum treatment can be performed properly 
throughout and leads to a finite unitary matrix. 

Multichannel quantum defect theory (MQDT) [11-15] and its classical limit 
[16-18] will be shown to provide the framework for such a paradigm. Indeed 
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we shall see that a simplified model of the Rydberg molecule allows to con- 
struct a classical Poincare map on the unit sphere, whose exact quantization 
is provided by MQDT. Thus the result is necessarily entirely quantal, exactly 
unitary and for finite h given in terms of a finite matrix. We shall show that 
the results commonly derived for MQDT are directly properties of the unitary 
representation of this classical map as obtained by MQDT. 

After a short description of the model for a Rydberg molecule and the simpli- 
fication introduced in Ref. [16], we proceed to give the quantum map for this 
case explicitly. We illustrate the two important aspects of our result by two 
applications. First the new interpretation allows modifications of the MQDT 
method, that prove particularly effective in near integrable systems. Second 
we proceed to show by way of examples that the properties of this map are 
relevant to the study of chaos and order in this system. 

Simplifying to the most basic case, these molecules can be viewed as a ro- 
tating system with positive charge and cylindrical symmetry that binds one 
electron in an orbit that is at large distances hydrogenic. The classical limit 
of the MQDT is the following classical model [16]: The motion is composed of 
two consecutive steps, (i) when the electron is far from the molecular core (i.e. 
most of the time for a Rydberg electron) it feels only the Coulomb part (— 1/r) 
of the potential. Its orbit is hydrogenic and its angular momentum L is fixed 
in the laboratory reference frame. Meanwhile the core rotates freely with an 
angular momentum N which is also fixed in the laboratory frame. The total 
angular momentum J = L + N is always conserved. In the molecular reference 
frame, the OZ axis is the cylindrical symmetry axis of the core. The core 
angular momentum N points in a perpendicular direction, taken as the OX 
axis. The angles 6 e and <p e are the polar and azimuthal angles respectively of 
the electronic angular momentum L in this frame. During this step, L rotates 
freely around the OX axis, (ii) during the so called "collision" step, the elec- 
tron senses also the cylindrically symmetric short range part of the potential 
of the core. Aside from the energy and J, the projection of L onto the core 
axis A = L cos 9 e is conserved due to the cylindrical symmetry of the core. We 
will add an extra, simplifying, hypothesis, namely that the magnitude L of L 
remains constant [16]. This is justified for Rydberg Molecules at least for small 
L's, but the classical and quantum map with this approximation exist for all 
L. Thus the collision can be described by a # e -dependent rotation of L around 
the core axis. The simplest form of this rotation compatible with the symme- 
try is [16]: 5<f e = K cos^e, where K is a coupling constant. This simplification 
is not essential. Notice further that the conservation of the total angular mo- 
mentum J implies that the molecular core feels a simultaneous recoil which 
changes the direction and magnitude of N. This change of N in turn entails 
a change of the rotational energy En of the core and because of conservation 
of total energy a change of the energy E e of the electron. This exchange of 
energy makes this model much richer than the kicked spin model [19] (which 
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is its limit when L <^ J, where this recoil can be neglected). In particular 
the energy of the electron may become positive after the collision, allowing 
to treat on equal footing bound and unbound (ionized) states. Possibility of 
chaotic motion comes from the conflict of these two steps, which consist of 
two rotations around distinct axes with different laws. The classically chaotic 
case can be obtained by increasing the coupling K. Near-integrable cases can 
be obtained for small coupling or at resonance, i.e. when the period of the 
electron is a multiple of half the period of rotation of the core. 

The quantum problem is solved by using the MQDT. The configuration space 
is divided by a sphere of radius r in a collision (r < r ) and an asymptotic 
region (r > r ) for the motion of the electron. r is of the order of the core 
size and is chosen such that in the asymptotic region the potential acting 
on the electron is only Coulombic, whereas in the collision region it feels 
both Coulomb and cylindrical potential. The wave functions for both regions 
are joined appropriately at r = r . The conflict between the two motions is 
expressed in quantum mechanics by the existence for the wave function of two 
bases with different good quantum numbers (in addition to J, J Z1 L and total 
energy E). At short distance the Born Oppenheimer basis is appropriate. 
The Rydberg electron is strongly bound to the core, thus quantized in the 
molecular reference frame and the additional good quantum number is A. At 
long distances the collision basis is appropriate. Here the electron is uncoupled 
from the core and the angular momentum N of the core remains a good 
quantum number. The collision is described by phase shifts fj,\, which are 
identical to collision phase shifts if the total energy is positive enough for all 
channels to be open, and which are related to the classical K parameter by 
//a = IM> ~ (K/AttL)A 2 . 

We focus on the completely bound situation, when total energy is low enough 
for the electron energy to be always negative whatever the value of N within 
the allowed range [J — L, J + L]. Demanding that the electron wave function 
goes to zero when r — > oo leads to demanding that the following determinant 
vanishes [15], i.e. 

det S = det {Una sin(7r(/ZA + vn))} = 0, (1) 

where the unitary U matrix given by 

U N , A =< L, -A, J, A\N, > (-1)^+ A (2 - 5 A , 0) 1 / 2 . (2) 

relates the two conflicting bases through a Clebsh Gordan coefficient. The 
principal (non integer) action vn{E) of the Coulomb electron, is related to the 
electron energy through E = E N + E e = N(N + 1)/ (2/) — 1/ (2v N ) , where I is 
the moment of inertia of the core, and we use atomic units (e = m = h = 1). 
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Corresponding wave functions are the eigenkets of S for the eigenvalue zero, 
labeled by A in the Born Oppenheimer basis: 

S\A A ) = 0. (3) 

Similarly the eigenbras (B N \ of S (eigenkets of S*), labeled by N describe the 
corresponding wave function in the collision basis. 

To proceed notice first that S is the imaginary part of a complex unitary 
matrix 

E = C + iS = {U NA exp(iir(nA + v N ))} . (4) 

This non symmetrical matrix maps the iV basis on the A basis: it is half 
the QPM we look for, and describes the motion between apogee and perigee. 
From it we can construct two QPM, which, by construction, turn out to be 
symmetrical unitary complex matrices. 

E*E operates in Born Oppenheimer A space, and is the exact quantization 
of a classical Poincare map. The latter is nearly the classical map used in 
refs. [16,17]. That map on the unit sphere described the position of L in 
the molecular frame immediately after the collision, whereas the present map 
describes the position of L in the middle of the collision (perigee). That this 
matrix is the T matrix defined by Bogomolnyfl] to quantize a Poincare map 
will be shown by proving that the eigenvalues and eigenfunctions of the entire 
system result from 

det(l-T(£ n )) = 0, (5) 

i.e. Bogomolny's equation for the quantized energy E n . Indeed at quantized 
energies given by Eq. (1) 

E*E|A A ) = 1|A A ), (6) 

which implies (5). To prove this key point first notice that unitarity of E, 
namely E f E = (C* - iS*)(C + iS) = I, leads to C*C + S*S = I and C'S = S*C. 
Then E*E = (C* + iS*)(C + iS) — I — 2S*S + 2iC*S, so that if S\A A ) = then 
BE\A A )=I\A A ). 

Conversely EE* operates in collision N space and corresponds to a classical 
map between apogee and apogee in the laboratory frame. 

We will now compare the traditional way of solving MQDT to the one implied 
by our QPM. The traditional way is to look for zeros of the determinant of the 
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non-symmetric matrix S of Eq. (1), all of whose elements depend in a complex 
way on energy through v^{E). It is computed by a LU or a SVD method 
followed by a root searching algorithm [21,16]. The present method is to look 
for eigenphases of E*E or EE*. They are computed by the diagonalization of 
a symmetrical unitary matrix. This is efficient and unproblematic because 
it diagonalizes in an orthonormal basis. Finally we search the zeros of the 
resulting eigenphases. 




251.0 251.2 251.4 251.6 251.8 252.0 
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Fig. 1. Comparison of Determinant and Eigenphases methods to solve the MQDT 
problem. Abscissa axis is energy on a scale vn{E) for the middle N = J value. 

The situation is sketched in Fig. 1. The search for zeros of eigenphases which 
vary nearly linearly with energy is obviously much simpler than the search of 
zeros of a determinant which is sometimes nearly tangent to the horizontal 
axis. Moreover, in this case of near tangency we had problems to converge 
the wave functions because the eigenfunction switches between two nearly 
orthogonal values in a very narrow energy range [17], requiring the use of the 
more efficient but slower SVD algorithm. Such a situation occurs frequently in 
nearly integrable cases, due to the lack of level repulsion. The diagonalization, 
on the contrary, gives always correctly the orthogonal eigenfunctions even if 
eigenvalues are very close. 

The implications of this procedure for the study of quantum chaos is of great 
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Fig. 2. Eigenphases in a chaotic case. 

interest. Figure 2 displays a chaotic situation. Comparing eigenphases in the 
near integrable (Fig. 1) and chaotic (Fig. 2) situations, we see that in the first 
case the phases as a function of energy display avoided crossings of straight 
lines running at different angles while for the second case these lines are prac- 
tically parallel. This shows that the drift of eigenphases as a function of energy 
displays presence or absence of level repulsion or spectral rigidity more obvi- 
ously than the energy levels themselves. This consideration is important in 
relation with the theory put forward by two of us [22], which relates Random 
Matrix Properties of eigenvalues of a quantum system to properties of invari- 
ance under canonical transformations of the structure of the corresponding 
classical system (structural invariance). This theory was developed for maps 
such as the scattering map, the stroboscopic map or the Poincare map leading 
to results about their unitary representations, i.e. about eigenphases. To trans- 
fer statistical properties of eigenphases of the QPM to energy eigenvalues, it is 
necessary that the drift of eigenphases as a function of energy be nearly paral- 
lel for all phases. Analytic evidence that this must be true for chaotic systems 
is given in [22,23], in agreement with the numerical results shown in Fig. 2. To 
put this on a more quantitative basis we display on Fig. 3 histograms of the 
distributions of the velocities of the eigenphases curves. The chaotic case has 
a narrow distribution while the integrable one shows long tails. The present 
study is much more convincing in that respect than another one for billiards 
using directly Bogomolny theory [23], which is garbled by the non unitarity of 
the T matrix for finite h. 

The difference seen in Fig. 3 for the near integrable and the chaotic case is 
remarkable and it is thus tempting to consider the eigenphase velocity distri- 
bution as another signature of classical chaos in quantum mechanics. While 
the narrow distribution for chaotic systems is typical, the tails for the inte- 
grable systems are not generic: they depend on the details of the partition of 
phase space by separatrices. In other examples [23] tails have different shapes. 
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Fig. 3. Histograms of the velocities of the eigenphases curves for a given energy. 
Top: chaotic case, bottom: integrable case. These histograms correspond to the 
same classical parameters as in Figs. 1 and 2, but H has been divided by 25 to 
increase statistics 

Summarizing: in this paper we proposed to interpret the multichannel quan- 
tum defect theory as a quantum Poincare map. This was implemented in detail 
in the approximation that the absolute value of the electron angular momen- 
tum is conserved, but it is quite clear that it is true in general. The new 
interpretation allows for a more stable and efficient way to find solutions to 
MQDT, particularly for near-degenerate levels. Beyond that, the approximate 
system has been used to exemplify quantum features of classically chaotic 
systems. This can be now extended to use MQDT as a paradigm for QPM. 
Indeed a study of the velocity distribution of eigenphases confirms properties 
expected from other studies. 

This work was partially supported by DGAPA (UN AM) project IN... and the 
CONACYT grant 25192-E 
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